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We introduce a method for finding general solutions of third-order nonlinear dif- 
ferential equations by extending the modified Prelle-Singer method. We describe 
a procedure to deduce all the integrals of motion associated with the given equa- 
tion so that the general solution follows straightforwardly from these integrals. The 
method is illustrated with several examples. Further, we propose a powerful method 
of identifying linearizing transformations. The proposed method not only unifies all 
the known linearizing transformations systematically but also introduces a new and 
generalized linearizing transformation (GLT). In addition to the above, we provide 
an algorithm to invert the nonlocal linearizing transformation. Through this proce- 
dure the general solution for the original nonlinear equation can be obtained from 
the solution of the linear ordinary differential equation. 
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1. Introduction 

In a previous paper (Chandrasekar et al. 2005) we have discussed the complete inte- 
grability aspects of a class of second-order nonlinear ordinary differential equations 
(ODEs) through a nontrivial extension of the so called Prelle-Singer (PS) (Prelle & 
Singer 1983; Duarte et al. 2001) procedure. We have illustrated the procedure with 
several physically interesting nonlinear oscillator examples. We have also developed 
a straightforward algorithmic way to transform the given second order nonlinear 
ODE to a linear free particle equation, if it is linearizable. 

One of the questions raised at the final stage of our earlier work (Chandrasekar 
et al. 2005) was that what are the implications of the novel features which we 
introduced in the extended PrcUe-Singer procedure to obtain the second constant 
of motion (in the case of second order ODEs) to third and higher order ODEs. To 
have a closer look at the problem let us recall our earlier work briefly here. We 
considered a second order ODE of the form ^ = g|*'|^'|!;| , P,Q & C[t,x,x], and 
explored two pairs of independent functions, say, Ri and Si, i = 1,2, associated 
with the underlying ODE. These functions are nothing but the integrating factors 
and null forms, respectively. Once these two pairs of functions are determined (by 
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solving an overdetermined system of first order PDEs) then eacli pair leads to 
an independent integral of motion, which can then be used to find the general 
solution for the given equation. Thus instead of integrating the first integral and 
obtaining the general solution which is conventionally followed in the literature we 
implemented some novel ideas in the PS method such that one can construct the 
general solution for the given equation in a self contained way and, in fact, our 
procedure works for a class of problems. 

In the case of third order ODEs, one should have three independent integrals 
of motion in order to establish the complete intcgrability. To deduce these three 
integrals one should have three pairs of independent functions {Ri, Ui and 6*1), 
i = 1,2,3. When we extend the PS procedure to third order ODEs we find that 
the determining equations for the integrating factors and null forms provide cither 
one or two integrals of motion only straightforwardly. Again the hidden form of the 
functions (i?3, 1/3,83) should be explored in order establish the complete intcgrabil- 
ity of the given equation within the frame work of PS procedure. In this paper we 
describe a procedure to capture the required set of functions. With the completion 
of this task we formulate a simple, straightforward and powerful method to solve a 
wide class of third order ODEs of contemporary literature. 

We stress at this point that the application of PS procedure to third order 
ODEs is not a straightforward extension of the second order case. In fact one has 
to overcome many faceted problems. The first and foremost one is how to solve the 
determining equations in such a way that one could obtain three sets of independent 
functions, namely, {Ri, Ui, Si), i = 1,2, 3, in a systematic way. In the present case 
we have six equations for three unknown functions (in the case of second order 
equations we have three equations for two unknowns). We overcome this problem 
by adopting suitable methodologies, the details which we present in §3. Another 
obstacle one could face in higher order ODEs, at least in some cases, is that one 
may be able to get only one integral of motion and in this situation how one would 
be able to generate the remaining integrals of motion from the first integral is also 
tackled by us in this paper. 

Our main goal, besides the above, is to bring out a novel and straightforward 
way to construct linearizing transformations for third order ODEs. The latter can 
be used to transform the given third order nonlinear ODEs to a linear equation. We 
note that unlike the second order equations, the third order ODEs can be linearized 
through different kinds of transformations, namely, invertible point transformation 
(Ibragimov & Meleshko 2005), contact transformation (Bocharov et al. 1993; Ibrag- 
imov & Meleshko 2005) and generalised Sundman transformation (Berkovich & 
Orlova 2000; Euler et al. 2003; Euler & Euler 2004). 

In this paper we introduce a new kind of transformation which can be effectively 
used to linearize a class of nonlinear third order ODEs. In fact, one can linearize 
certain equations only through this transformation alone and not by the known 
ones in the literature. We call this transformation as generalized linearizing trans- 
formation (GLT). We note that generalised Sundman transformation is a special 
case of this transformation. In the generalised Sundman transformation the new 
independent variable is a nonlocal one and so eventhough one is able to transform 
the given nonlinear third order ODE to a linear one due to the nature of nonlocal 
independent variable it is not easy to write down the general solution. In the case 
of generalized linearizing transformation both the new dependent and independent 
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variables contain derivative terms also in addition to the independent variable being 
nonlocal. Even for this general case, in this paper, we succeed to present an efficient 
algorithm to deduce the general solution. 

Another fundamental problem regarding linearization is how to deduce the 
linearizing transformations systematically. Generally, Lie symmetry analysis and 
direct methods arc often used to deduce the point and contact transformations 
(Steeb 1993; Olver 1995; Bocharov et al. 1993; Bluman & Anco 2002; Ibragimov 
& Meleshko 2005). In this work we propose a simple and straightforward method 
to deduce linearizing transformations and we derive them from the first integral. 
Our method of deducing linearizing transformations has several salient features. Ir- 
respective of the form of the linearizing transformation (point/contact/generalised 
Sundman transformation), it can be derived from the first integral itself. We also 
note that one can also linearize a third order ODE to the second order free particle 
equation through our method. An added advantage of our method is that suppose 
a given equation is lincarizablc through one or more kinds of transformations then 
our procedure provides all these transformations in a straightforward way and as far 
as our knowledge goes no such single method has been formulated in the literature. 

The plan of the paper is follows. In the following section wc extend the PS proce- 
dure to third order ODEs and indicate new features in finding the three independent 
integrals of motion. In §3, we describe the methods of solving the determining equa- 
tions and how one can obtain compatible solutions from them. Wc illustrate the 
procedure with several examples. In §4, we propose a powerful method of identi- 
fying linearizing transformations. This method not only brings out all the known 
transformations systematically but also a new generalized linearizing transforma- 
tion (GLT) for the third order ODEs. We emphasize the validity of the method with 
several illustrative examples arising in different areas of mathematics and physics. 
We present our conclusions in §5. 



2. Prelle-Singer method for third order ODEs 

(a) General Theory 
Let us consider a class of third order ODEs of the form 

x = ^, P,QGC[t,x,x,x], (■=^) (2.1) 

where over dot denotes differentiation with respect to time and P and Q arc poly- 
nomials in t, X, X and x with coefficients in the field of complex numbers, C. Let 
us assume that the third order ODE (2.1) admits a first integral I{t,x,x,x) = C, 
with C being constant on the solutions, so that the total differential of I gives 

dl = Itdt + Ixdx + hdx + hdx = 0, (2.2) 

where each subscript denotes partial differentiation with respect to that variable. 
Equation (2.1) can be rewritten as ^dt — dx = 0. Now adding the null terms 
U{t,x,x,x)xdt—U{t,x,x,x)dx and S{t,x,x, x)xdt—S{t, x, x, x)dx to this we obtain 
that, on the solutions, the 1-form 

P 

{—+Sx + Ux)dt-Sdx-Udx-dx = 0. (2.3) 
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Looking at equations (2.2) and (2.3) one can conclude that, on the solutions, 
these two forms are proportional and the form of equation (2.3) is equivalent to 
equation (2.2) except for an overall multiplication factor. Thus multiplying equa- 
tion (2.3) by the factor R{t,x,x,x) which acts as the integrating factor for (2.3), 
we have on the solutions that 



dl = R{(f) + Sx + Ux)dt - RSdx - RUdx - Rdx = 0, 



(2.4) 



where (j) = P/Q. Comparing equations (2.2) with (2.4) we have the following rela- 
tions, on the solutions. 



It = R{(l) + Sx + Ux), 4 



-RS, 



h = -RU, 



-R. 



(2.5) 



Now imposing the compatibility conditions, 7^ = Ixt, hi = ht, hx = ht, hx = 
Ixx, Ixx = Ixx, Ixx = Ixx, which exist between the equations (2.5), wc have the 
following equations which constitute three determining equations ((2.6)-(2.8) given 
below) for the functions S, U and R along with three constraints ((2.9)-(2.11) given 
below) that they need to satisfy. 



D[S] = -(t>x + Scj>i + US, 



D[U] 

D[R] : 
Rx 

RxS ■ 

Ri ■ 



-(bx + U<t>x-S + U' 
~R{U + (l)x), 
RxS + RSx, 

— RSx + RxU + RUx; 

RxU + RUx, 



where 



D = 



d_ 

dt 



dx 



dx 



dx 



(2.6) 

(2.7) 
(2.8) 
(2.9) 
(2.10) 
(2.11) 

(2.12) 



The task of solving equations (2.6)-(2.11) can be accomplished in the following 
way. Substituting the given expression of (/) into (2.6)-(2.7) and solving them one 
can obtain expressions for 5* and U. With the known U, equation (2.8) becomes the 
determining equation for the function R. Solving the latter one can get an explicit 
form for R. Compatible solutions to equations (2.6)-(2.8) can also be obtained in 
alternate ways, the details of which are given in §3. 

Now the functions R, U and 5* have to satisfy an extra set of constraints, that 
is, equations (2.9)-(2.11). Suppose a compatible solution satisfying all the equa- 
tions have been found then the functions R, U and S fix the differential invariant 
I{t, X, X, x) by the relation 



I{t,x,x,x) =ri-r2- j |i?C/ + ^ I*^! ~ '^2]| 



dx 



R + 



dx 



ri - r2 - 



iif/ + — [ri - ra] } dx 



dx, (2.13) 



where 



R{(p + Sx + Ux)dt, r2 



RS 



dx 



ri ] dx. 
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Equation (2.13) can be derived straightforwardly by integrating the equations (2.5). 
Here it is to be noted that for every independent set {S,U,R), equation (2.13) 
defines an integral. 

(b) Exploring the eomplete form of R: Theory 

From the above discussion, it is clear that equation (2.1) may be considered 
as completely integrable once we obtain three independent sets of the solutions 
{Si, Ui, Ri), i = 1,2,3, which provide three independent integrals of motion through 
the relation (2.13). Here we note that since we are solving equations (2.6)-(2.8) first 
and then check the compatibility of this solution with equations (2.9)-(2.11), one 
often meets the situation that all the solutions which satisfy equations (2.6)-(2.8) 
need not satisfy the constraints (2.9)-(2.11) since equations (2.6)-(2.11) constitute 
an overdetermined system for the unknowns R, S and U. In fact for a class of prob- 
lems one often gets one or two sets of S, U, R which satisfy all the equations (2.6)- 
(2.11) and another set(s) {S, U, R) which satisfies only the first three equations and 
not the other, namely, (2.9)-(2.11). In this situation we find an interesting fact that 
one can use the integral derived from the set(s) which satisfies all the six equa- 
tions (2.6)-(2.11) and deduce the other compatible solution(s) {S,U,R) (definition 
of R follows). For example, let the set (^3, C/3, R3) be a solution of the determining 
equations (2.6)-(2.8) and not of the constraints (2.9)-(2.11). After analysing several 
examples we find that one can make the set (^s, U3, R3) compatible by modifying 
the form of Rs as 

R3 ^ F{t,x,x,x)R3, (2.14) 
where R3 satisfies equation (2.8), so that we have 

{Ft + xF^ + xFi + (j)Fi)R3 + FD[R3\ = -Fi?3(C/3 + (pi). (2.15) 

Further, if F is a constant of motion (or a function of it), then the first term on 

the left hand side vanishes and one gets the same equation (2.8) for R^ provided 
F is non-zero. That is, whenever is a constant of motion or a function of it then 
the solution to (2.8) may provide only a factor of the complete solution Rs without 
the factor F in equations (2.14). This general form of R3 along with S3 and U3 can 
now provide complete solution to equations (2.6) - (2.11) as discussed below. 

(c) Exploring the complete form of R: Method 

Now if the sets {Si,Ui,Ri), i = 1,2, are found to satisfy the equations (2.6) 
- (2.11) and the third set {S3,U3, R3) does not satisfy equations (2.9) - (2.11) 
then F may be a function of the integrals Jj, i = 1,2, derived from the sets 
{Si,Ui,Ri), i = 1,2. We need to find the explicit form of F{I\,l2) in order to 
obtain the compatible solution {SsjUsjRs). To do so let us find the derivatives of 
R3 with respect to x, x and x: 

Rsx = {F[hx + F2l2x)R3 + FR3X, Rsx = {Fihx + F2^2x)R3 + FR3X, 

Rzx = {F[hx + F!,hx)Rz + FR3X, (2.16) 
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where Fl = ^ and = Substituting equation (2.16) into equations (2.9)- 
(2.11), we have 

ihFj + f2F^) ^ F ^ {f,Fi+f,F^) ^ F ^ ifrFj + fsF^) ^ F ^ ^^.17^) 



where 



fl = {hx — hxSs), = {hx — hxSs), /s = {S3R3X + R3S3X — Rsx) 
fi = {hx — hxUa), fh = [hx — hxUs), /e = (f/s-Rsx + RaU^x — Rsx) 
fl = {Sahx — hxUs), fs = {S^hx — hxUa), 

/g = {R3U3X + U3R3X — R3S3X — SsR^x)- (2.17b) 



Equation (2.17a) represents an ovcrdctcrmincd system of equations for the un- 

dF ' 
ail ' 



known F . A simple way to solve this equation is to uncouple it for F[ (= ^) and 



F^ (= ^) and solve the resultant equations. For example, eliminating F'^ from 
equation ^2.17a) we obtain equations for F[ in the form 

RiF[ ^ ihfr, - J2f6) ^ ifsfs - f-ifo) ^ ifefs - ./:-../•») . . 

F (/1/5 - /2/4) (/1/8 - /2/7) (/4/8 - 75/7) " ^ ' ' 

On the other hand eliminating F{ from equation (2.17a) we arrive equations for F2 
in the form 

^ (/3/4 - fife) ^ ihfl - /1./9) ^ (/6./7 - /4/9) . . 

F (/2/4 - /l.f5) (/2/7 - ,/l./8) (/5,/7 - ^./s) ' ^ ' ' 

It can be easily cheeked that the compatibility of the right three expressions 
in equations (2.18) or (2.19) give rise to relations which are effectively nothing 
but the constraint equations (2.9)-(2.11) and so no new constraint is added now. 
Consequently, equations (2.18) and (2.19) can be written as 

|£=5(7i,72)i^ and ^ = h{h,l2)F, (2.20) 

respectively, where g{h,h) = and /i(Ji, /a) = Ws ulu'JflfJ, - ^^"^ 

can solve equations (2.20) and obtain the form of F{Ii,l2). This is demonstrated for 
several examples in the following sections explicitly. Once F is known wc can obtain 
the complete solution i?3 from which, along with ^3 and t/3, the third integral h 
can be constructed. Thus with the explicit forms of the three integrals of motion, 
the complete integrability of equation (2.1) is guaranteed. 

Finally, if the set {Si,Ui,Ri) alone is found to satisfy the equations (2.6) - 
(2.11) and the second set {S2, UijRi) also does not satisfy equations (2.9) - (2.11), 
then F may be a function of the integral Ii alone which was derived from the 
set (^i, Ui,Ri). We need to find the explicit form of F{Ii) in order to obtain the 
compatible solution {82,1/2, Ri)- To do so let us find the derivatives of R2 with 
respect to x, x and x: 

Rlx = F[lixR2 + FR2X, R2x = F[lixR2 + FR2X, R2x = F[lirj.R2 + FR2X, 

(2.21) 
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where = Substituting equation (2.21) into equations (2.9)-(2.11), we have 

^ = ^ = ^ = ^, (2.22a) 

F Si S3 S5 

where 

Sl = {hx — hxS2), S2 = {S2R2X + R2S2X — R2x), 

S3 = {hi - hxU2), S4 = {U2R2X + R2U2X - R2x), 

S5 ={S2hx-hxU2), S6 = {R2U2x + U2R2x-R2S2x-S2R2x). (2.22b) 

One can again check that the compatibiHty of the right three expressions in equation 
(2.22a) leads to a condition, which is deducible from (2.9)-(2.11), and so no new 
condition is introduced in reaUty. Then rewriting equation (2.22a) as 



dh ~ \R28j 



F = g{h)F, (2.23) 



and solving it, one can obtain an explicit expression for F. Once F is known we can 
construct the complete form of R2 from which, along with S2 and U2, the second 
integral of motion can be obtained. Once h and I2 are known, we can proceed to 
find the third compatible set (^3, Us, R3) as before and obtain the third integral /3 
to establish complete integrability. 

3. Methods of finding the expHcit form of R 

In §2 we outlined the method of solving the determining equations. However, in 

practice, it is difficult to solve equations (2.6)-(2.8) straightforwardly as they con- 
stitute a set of coupled first order nonlinear partial differential equations (PDEs). 
So one has to look for some intuitive ideas to solve these equations. We solve them 
and obtain the forms R, U and S in the following way. For this purpose, we observe 
the important fact that when we rewrite the coupled equations (2.6)-(2.8) into an 
equation for a single variable, namely, R, the resultant equation turns out to be a 
linear PDE. Then wc solve this "i? equation" with a suitable ansatz (say polyno- 
mial or rational in x). Once 'R' is found the remaining functions U and S can be 
easily deduced. 

As noted, rewriting equations (2.6)-(2.8) we arrive at a third order linear partial 
differential equation for R in the form 

D^[R]-D^[R4,x]+D[Rcf>x]-<f>xR = 0, = + + + '^g^' ^^'^^ 

Solving equation (3.1) with a suitable ansatz in x is relatively easier in many cases 
than solving equations (2.6)-(2.8). Once the explicit form of R is obtained, U can 
be deduced from equation (2.8) as 

C/ = -(^+<^x), (3.2) 

from which S can be fixed by using equation (2.7). Now if this set (5, U, R) forms 
a compatible set for the remaining equations (2.9)- (2.11) then the corresponding 
integral / can be found using equation (2.13). To illustrate this idea let us look into 
the following examples. 
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(i) Example:! 

Let us begin with a simple example, namely a linear third order ODE, 

X + Ax = 0, (3.3) 

where A is an arbitrary parameter. Substituting (p = — Aa; into (3.1), we get the 
following equation for R 

D^[R] + XR = o, ^ = | + .^|: + 4-^4- ^'-'^ 

We now assume an ansatz for R in the form 

R = a{t,x,x) + b{t,x,x)x, (3.5) 

where a and b are arbitrary functions of t, x and x. Substituting (3.5) into (3.4) 
and equating the coefficients of different powers of x to zero we get a set of linear 
partial differential equations for the variables a and h. Solving them one obtains 
three nontrivial solutions, 

i?i=-e^*, R2 = -{2x + kx-k'^x)e-^\ R3 = ~^{x + kx)e-''\ (3.6) 

where = A. Now substituting the form of Ri's ,i = 1,2,3, separately into the 
equation (3.2), we get 

{2k^x + kx + k^x) {x - ex) 

{2x + kx-k^x) ' {x + kx)- ^^■^> 

Now substituting the forms of C/j's, i = 1,2, 3, into (2.7) one can fix the forms of 
Si S) i —— Ij 2j 3, cis 

S^ = k\ - + ^-"^ ^ 53=_;tli±M. (3.8) 

' 2x + kx — k'^x ' {x + kx) 

As a consequence now we have three sets of independent solutions for the equation 
(2.6)-(2.8). Now we check the compatibility of these solutions with the remaining 
equations (2.9)-(2.11). 

We find that the solutions (^i, Ux,Ri) and (5*2, ^2) satisfy the equations (2.9) 
- (2.11), and become compatible solutions. Substituting the forms (5*1, Ui.Ri) and 
{S2, U2, R2) separately into equation (2.13) and evaluating the integrals we get 

I-^ = ^(^-kx + k^x^ e*^*, (3.9) 

72 = ^ + fc^j-a _^ ^4^2 ^ _ ^2^j _^ k^xx^ e^''^. (3.10) 

However, the set {83,113, R3) docs not satisfy the extra constraints (2.9)-(2.11) 
which means that the form of R3 may not be the 'complete form' but might be a 
factor of the complete form. As mentioned in §2, in order to recover the complete 
form R3 one may assume that R3 = F(7i,/2)i?3, where F{Ii,l2) is a function of 
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the integrals Ji and Now we have to determine the form of F{Ii,l2) expUcitly 
and for this purpose we proceed as follows. Substituting 



Rs = F{h,h)Rs 



F{h,h), 



(3.11) 



into equations (2.18) and (2.19), along with (2.17b), we obtain two equations for F 
as, 

F[ = 0, hF!^ + F^0, (3.12) 

where F{ and F2 denote partial derivatives of F with respect to Ii and I2, respec- 
tively. Upon integrating (3.12) we get, F = j^, (the integration constants are set 

to zero for simplicity) which fixes the form of Rs as 



i?3 



R3 



{x + kx) 



2 [x^ + k'^x^ + fc^a;^ + kxx — k'^xx + k^xx) ' 



(3.13) 



Now one can easily check that this set (^3, U3, R^) is a compatible solution for the 
equations (2.6)-(2.11) which in turn provides /s through the relation (2.13) as 



h - - 



t + tan 



X — kx ~ 2k X 
\/3{x + kx) 



(3.14) 



Using the explicit forms of the integrals h, I2 and I3, the solution to equation (3.3) 
can be deduced directly as 



x{t) = /le"*^* + ^/he^* cos 



/\/3fc 



t + F 



(3.15) 



The result exactly coincides with the solution presented in Polyanin & Zaitsev 
(2000). 



(ii) Example: 2 

The applicability of this method to nonlinear ODEs can be illustrated by con- 
sidering an equation of the form 

"2 • " 

X = — + —. (3.16) 

X X 

Equation (3.16) is a sub-case of the general form of a scalar third order ODE 

which is invariant under the generators of time translation and rcscaling (Feix et 
al. 1997; Polyanin & Zaitsev 2000). A sub-case of equation (3.16), namely, 'x — 
= 0, has been considered by both Bocharov et al. (1993) and Ibragimov & 
Mclcshko (2005) to show that it can be linearized to a linear third order ODE 
through a contact transformation. On the other hand, Euler & Euler (2004) have 
considered the equation x — ^ = and showed that it can be linearized through 
Sundman transformation (sec §5c below). Here wc consider the combined form 
(3.16) and derive integrating factors, integrals of motion and the general solution 
of this equation. Further, we show that the equation (3.16) itself can be linearized 
by the generalized linearizing transformation (GLT) (see §5d below). 
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As before, substituting ^ = ^ + ^ into (3.1), we get the following linear partial 
differential equation for R as 

D^R] - Z)^[(| + l)R] - D[i^ - l)R] + §R=0, (3.17) 

where D is defined by equation (2.12). Now substituting the ansatz (3.5) into (3.17) 
and proceeding as before, we get 

R. = -^, = ij3 = !^!z_E(i±i^. (3.18) 

XX X 2xx'^ 

Following the ideas given in Example 1, one can deduce the corresponding forms of 
SiS and J/j's, i = 1,2, 3, as 

51 = --, U, = ~, (3.19) 

X X 

52 = -, = — - ?, (3.20) 

X XX 

53 = a:a;(a; + tx) ^ xx{2x + tx) 

x{—tx^ + x{x + tx))^ x{tx'^ — x{x + tx))' 

The solutions {Si,Ui,Ri), i = 1,2, satisfy the constraints (2.9)-(2.11) so that they 
lead to first and second integrals of the form 

h = —, l2 = : . 3.22 

XX X 

In the present case also the set {8-^,11^, R^) docs not satify the extra constraints 
and so one has to explore the complete form of R^- To do so we proceed as before 
and obtain the forms of F and R^ bs, F = ,} , and 

V Jl J2 

2(^Iil2)xx^ 

where the explicit forms of Ii and I2 are given in equation (3.22). Now one can check 
that the set {S3, U3, R3) satisfies all the six equations (2.6)-(2.11) and furnishes the 
third integral I3 through the relation (2.13) as 



h = -l-{Vhh)t + tan-\^x. (3.24) 
^ \ ^2 



Using the explicit forms of the integrals Ii, I2 and I3, the solution to equation (3.16) 
can be deduced directly as 



x{t) = \ -f- tan 
V h 



^(Viy2i + 2/3) 



(3.25) 



As can be seen from equation (3.23) the complete compatible solution R3 has x 
term which appear inside the square root sign. This form of can be explored 
only by making a suitable ansatz. Moreover one may also face more difficulties in 
solving the determining equations (2.6)-(2.8). In such complicated situations the 
complete solution R can be obtained by using our procedure. 



Article submitted to Royal Society 



Integrability and Linearization 



11 



(iii) Example:3 

Let us consider a Chazy class of equation of the form (Halburd 1999; Mugan & 
Jrad 2002; Euler & Euler 2004) 

X + Aaxx + iax^ + Qa^x^x + a^x'^ = 0. (3.26) 

Substituting (p = —{Akxx + Skx"^ + Qk'^x^x + k^x'^) into (3.1), we get 

D^[R] + AaD\xR] - QaD[(x + ax'^)R] + 4q!(x + a^x^ + Zaxx)R = 0. (3.27) 

To solve the above equation wc assume the same ansatz (3.5) for the variable R. 
However, substituting this ansatz into equation (3.27) and solving the resultant 
equation we obtain only trivial solution. So we seek a rational form of ansatz for R 
in the form 

^ ^ a{t,x,x)+h{t,x,x)i 
{c{t, x, x) + d{t, X, x)xY ' 

where r is an arbitrary number. Substuting the equation (3.28) into equation (3.27) 
and solving the resultant partial differential equations we get 

^ X + ax^ ^ t{—2x + atx^ + tx) 

{a^x^ + 3axx + xY ' {a^x'^ + iaxx + i)^ ' 

_ ^(3 - 2,atx + aH'^x'^ + at^x) 
{a^x^ + Zaxx + xY 

The corresponding forms of SiS and J/j's i = 1,2, 3, are 

^a^x^ - X a{a'^x^ + 3x'^ - 2xx) 

^1 " — 2 I ■ ' ^ 2~r^ ' 

ax'^ + X ax'' + x 

_2x- Qatx^ + 2a^t^x^ ~ t^x 

- t{-2x + atx^+tx) ' ^^-^^^ 



52 = 

53 = 



2ax^ - AaHx^ + aH^x^ -2x + Sat^x^ + 2tx{l - atx) 

t{-2x + + tx) ' 
3 - 12ate + Qa^tl^x^ - 2a^i^x^ + at^x 
t(3 — 'iatx + a'^t'^x'^ + at'^x) ' 
(a(12ate2 - Qa^t^x'^ + a^t^x^ + ■it{2x + at^x^ + tx) - 2x{i + at^x))) 
i(3 - "iatx + aH'^x'^ + at'^x) ' 



Now we proceed to find the integrals of motion. First we note that the functions 

(5i, Ui, Ri) satisfy the constraints (2.9)-(2.11) and hence they are compatible. Thus 
substituting them into (2.13) the first integral I\ is fixed easily as 

a.T^ + X ,„ „^ . 

h = -t+ — ; -. (3.31 

a-^x"^ + iaxx + X 

However, the set (<S'2, ?72, -R2) (and so also {Sz,Uz,R3)) does not satisfy the extra 
constraints (2.9)-(2.11) which means the form of i?2 niay not be the 'complete form' 
but might be a factor of the complete form. As mentioned in §2, in order to recover 
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the complete form R2 one may assume that R2 = F{Ii)R2. Here F{Ii) is a function 
of the integral I\. Now we have to determine the form of F{Ii) explicitly and for 
this purpose we proceed as follows. Substituting 

4 = F(/i)ii2 = gi^^H^m), (3.32) 

into equation (2.22a), we obtain IiF[ + 2F ~ 0, where F[ denotes differentiation 
with respect to h. Upon integrating the latter we get F = /j"^, (the integration 
constant is set to zero) which fixes the form of i?2 as 

i?2 = ^(-2x + atx"^ + tx) 

{ax^ — a^tx^ + X — iatxx — tx)^ 

Now one can easily check that this set (52, U2, R2) is a compatible solution for the 
equations (2.6)-(2.11) which in turn provides I2 through the relation (2.13) as 

^ _ -2atx'^ + aH^x^ + x{2 + 2>o£-x) + t(-2a; + fx) 
^ (ax^ — (xHx^ -\- X — Zoixx — tx) 

As in the previous examples, the set (^3, U3, R3) does not satisfy the constraints 
(2.9)-(2.11) and hence one should seek a complete form for R3, which we denote as 
R3, in the form 

i?3 =F{h,ms=F{h,l2y-^^^^^^^^. (3.35) 

Substituting (3.35) into equations (2.18) and (2.19), we obtain the following equa- 
tions for F, that is, hFi + 2F = 0, F^ = 0, where again F{ =^ and F^ = 
Upon integrating the equations we get the explicit form of F as = which in 
turn fixes the form of i?3 as 

_ <(3 - Satx + kH'^x'^ + at^x) 
3a(ax^ — aHx^ + x — Satxx — tx)'^ 

Now the set {S3, U3, R3) satisfies all the six equations (2.6)-(2.11) and the relation 
(2.13) gives the form of third integral I3 as 

^ _ (6 + Sd-f-.r- - n-'/-'.]:-' + 3n/-i: - 3n/,.i:(2 + n/.-j-j - (\l:^x) 
^ 6a{ax'^ — aHx^ + x — 3atxx — tx) 

Thus we have obtained the explicit forms of the integrals Ii , I2 and 73 and hence 
the solution to equation (3.3) is obtained directly as 



^+ht + hh 
— -r all ^ + ahht + hh 



<i) - o,^ 'rV"^'^:\" . r. - (3-38) 



Interestingly one can derive the solution (3.38) through an alternate way also. 
For example, instead of solving the 'i? equation' with rational ansatz, one can look 
for equations in other variables, that is, either in U or in S. For example, from 
equation (2.7) we get 

S = -{D[U] U(i>i - f/2). (3.39) 
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Substituting equation (3.39) into the equation (2.6) we get a nonlinear PDE for U: 

D^[U] - 3UD[U] - UD[cPi] + D[cPi] - cj>, - (t>i<i>x + ^lU 

+2(/)aC/2 - (f>^U + = 0. (3.40) 

Now one can look for solutions of equation (3.40) with polynomial in x. Once 
U is known one can make use of equations (3.39) and (2.8) to get the forms of 
corresponding S and R respectively. It turns out that for some cases, like the present 
example (3.26) (sec the actual forms of U in equation (3.30)), solving equation (3.40) 
is easier than solving equation (3.1) to solve. However, this can be decided only by 
actual calculation. 



4. Linecirization 

In §2 and §3 we discussed the complete integrability of third order ODEs by inves- 
tigating sufficient number of integrals of motion. However, one can also establish 
the complete integrability of the given nonlinear ODE by transforming it into a 
linear free particle second order ODE or into a third order linear ODE of the form 
= w'" = 0. Unlike the second order ODEs, the third order nonlinear ODEs 
can be linearized through a wide class of transformations, namely, invertible point 
transformation (Steeb 1993; Ibragimov & Meleshko 2005), contact transformation 
(Duarte et al. 1994; Bocharov et al. 1993; Ibragimov & Meleshko 2005), generalised 
Sundman transformation (Berkovich & Orlova 2000; Euler et al. 2003; Euler & 
Euler 2004) and their generalizations. In the following, we describe a procedure to 
deduce these transformations from the first integral itself and illustrate our ideas 
with relevant examples. 

(a) Transformation from third order nonlinear ODEs to second order free particle 

equation 

Let us suppose that the ODE (2.1) admits a first integral 

h=F{t,x,x,x), (4.1) 



where F is a function o{t,x,x and x only. Now extending our earlier proposal for 
second order ODEs (Chandrasckar et al. 2005) to the third order equations (2.1), 
let us split the function F into a product of two functions such that one involves a 
perfect differentiable function ^Gi{t, x, x) and another function G2{t, x, x, x), that 
is, 

h=F f „ . ... ^G^{t, X, x)\ . (4.2) 

\G2{t,x,x,x) dt ) 

Suppose G2 is a total time derivative of another function, say, z, that is dz/dt = 
G2{t,x,x,x), then (4.2) can be further rewritten in the form 

/ 1 dGi\ fdGi\ , , 
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Now identifying the function Gi{t,x,x) = w as the new dependent variable, equa- 
tion (4.3) can be recast in the form 7i = F[^). In other words, we obtain 

= ^' ^'-'^ 

where /i is a constant, from which one can get = 0. Rewriting w and z in terms 
of old variables, namely, 

w = Gx{t,x,x), z= G2{t' ,x,x,x)dt' , (4-5) 



wc can get a linearizing transformation to transform the third order nonlinear ODE 
into the second order free particle equation. Then to deduce the general solution 
one has to carry out one more integration. 

(6) Transformation from third order nonlinear equation (2.1) to third order linear 

ODE w'" = 

Next, we aim to transform equation (2.1) into a third order linear ODE and so 

we try to rewrite the first integral as a perfect second order derivative. We note 
that this can be done when one is able to evaluate the following integral explicitly, 



w 



I Gi{t',x,x)G2{t',x,x,x)dt' = G3{t,x,x), (4.6) 

J o 



where Gi and G2 are defined as in equation (4.2). Note that in the function G3, 
the X dependence has been integrated out. The reason for making such a specific 
decomposition is that in this case equation (4.2) can be rewritten as a simple 
second order ODE for the variable w (see equation (4.8) below). We pursued a 
similar kind of approach in the integrable force-free Duffing van der Pol oscillator 
equation (Chandrasekar et al. 2004) which has now been generalized in the present 
case. Here one can rewrite (4.1) as a perfect second order derivative as follows. 
Differentiating (4.6) with respect to t, we obtain ^ = G1G2. Rewriting the left 
hand side in the form 

and using the identities already used in equation (4.5), namely ^ = G2 and Gi = w, 
in equation (4.7) the latter becomes 

dw , , „, 

— =w. (4.8) 
dz 

Differentiating (4.8) with repect to z and using the identity ^ = h (vide equation 
(4.4)), one gets 4-t = h- In other words, we have 



dz 

dz' 

d^w 



,,3-0, (4.9) 



so that w and z are the required transformation variables. 
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(c) The nature of transformations 

In the previous two subsections, §4 a and §4 6, we demonstrated how to con- 
struct Hnearizing transformations from the first integral and how they effectively 
change the given third order nonhncar ODE to cither second or third order linear 
equation. Depending upon the explicit form of the transformations we can call them 
as point, contact, generalised Sundman or generalized linearizing transformations 
(GLT). To demonstrate how different kinds of transformation arise let us consider 
the transformation, w = G3, z = J^G2dt' (vide equations (4.6) and (4.5)), which 
takes the given nonlinear ODE into a linear equation. Now, in the above transfor- 
mation, suppose z is a perfect diffcrontiable function and w and z do not contain 
the variable x and x, then we call the resultant transformation, namely, w = fi{x, t) 
and z = f2{x,t), as a point transformation. Further, if the transformation admits 
the variable x also explicitly then it becomes the contact transformation. In this 
case, we have w = fi{t,x,x) and z = f2{t,x,x). On the other hand, if ii) = G3{t,x) 
and z = J* G2{t' ,x)dt' then the transformation is said to be a generalised Sund- 
man transformation. Note that in the latter case the new independent variable is 
in an integral form. Besides the above, as pointed out earlier, we find that there 
exists another kind of transformation, namely, generalized linearizing transforma- 
tions (GLT) , in which the new dependent and independent variables take the form 
w = Gs{t,x,x) and z = J^G2{t',x,x,x)dt', respectively. 

(d) Connection between the functions Gi, G2 and G3 

Finally we explore the connection between the functions Gi, G2 and G3. As we 

have seen above, for the given equation to be linearizable it shoiild be transformable 
to the form (4.8). Rewriting the latter in terms of the variables, Gi,G2 and G3 we 
get, 

Gst + xG^x + iGa± = Gi(t,x,±)G2(t,x,x,x). (4.10) 

Note that the left hand side of equation (4.10) contains the variable x linearly. So 
the right hand side should also be linear in x. Consequently, we can write 

G2{t,x,x,x) = G2i{t,x,x)x + G22{t,x,x). (4-11) 

Using (4.11) we can rewrite equation (4.10) in the form 



Gat + XG3, 1 



XG3 



Gi. (4.12) 



Since the right hand side is independent of x, we have from (4.12) that 

G2i{G3t + xGsx) = G22G3x and Gi = (4.13) 

It may be noted that a similar condition has been derived by Bocharov et al. 
(1993) and Ibragimov & Meleshko (2005) for the case G2 is a perfect differential 
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function. In other words, our procedure indicates that more generahzed transfor- 
mations are possible in the case of third order ODEs. By imposing the condition 
(4.11) it becomes clear that whatever be the type of linearizing transformation, the 
new independent variable should be at the maximum linear in x. 



(e) Transformation to fourth and higher order linear ODEs 

In the above sub-sections §4 a , §4 & , §4 c and §4 d wc can concentrated only on 
transforming a nonlinear third order ODE either to a second order or third order 
linear equations only. However, one could also linearize certain third order nonlinear 
ODEs to fourth order linear ODE. For example the equation (3.26) is linearizable to 

,4 

a fourth order ODE of the form ^-r = 0, under the nonlocal transformation x = —. 

az^ ' aw 

This is not an isolated example and one can linearize a class equations through this 
procedure. Besides the above one can also consider linearizing transformations in 
which the new dependent variable, w, is a non-local one. This choice leads us to 
classify another a large class of equations which we leave it for future work. 



5. Application 

In this section we consider specific examples to demonstrate the method given in 
§4. 



(a) Example 1: Point transformation 

Let us consider a nontrivial example which was discussed by Steeb (1993) in the 
context of invertible point transformations, namely, 

a; H 3x h 2a; = 0. (5.1) 

XX 

The first integral, which can be obtained using the formulation in §2, can be written 
as 

^2 , " „ • \ -2t 



7i = l^a;^ +XX- xxj e"^\ (5.2) 

Rewriting (5.2) in the form (4.2) we get Ii — e~^-^{xxe~*), so that 

w = xie"*, z = e*. (5-3) 

As we noted earlier, one could transform (5.1) to the second order free particle 
equation, = 0, by utilizing the transformation (5.3). Integrating the equation 

= we get w = I\z + 12- Using (5.3) into this expression, the general solution 
of equation (5.1) can be obtained (after an integration) as 

x{t) = {he''' + he' + h)K (5.4) 

where h, i = li2,3, are the integration constants. 

Further, using the equation (4.6) we get w = J xxdt = Then we can directly 
check that the point transformation 

x'^ 

w=Y' z = e*, (5.5) 
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transforms the equation (5.1) to the form = 0. As we mentioned earher, since 
it) and z involve only x and t they are just point transformations. Integrating the 
linear equation = we get 

w = ^z^ + hz + l3. (5.6) 

Rewriting w"; and z in equation (5.6) in terms of the original variables x and t by 
using the transformation (5.5) we get the same solution as equation (5.4). 



(6) Example 2: Contact transformation 
Let us consider an equation of the form 

X = ^. (5.7) 

Bocharov et al. (1993) have shown that equation (5.7) can be linearized through 
contact transformation. However, the explicit linearizing transformation is yet to 
be reported which is also a difficult problem. Here we derive the explicit form of 
the linearizing transformation through our procedure. 

The first integral can be easily deduced using the results of §2 as 

• 2 

= (5.8) 



XX 



X d ,x. 



Rewriting (5.8) in the form (4.2) we get /i = — r(~) so that we have w = ^ and 

X at X ^ 
z = logx. The latter transforms equation (5.7) to the second order free particle 

,2 

equation — 0, so that w = Iiz + 12, where Ii and I2 are integration constants. 
Rewriting w and z in terms of the old variables we get x = (Ji log(i;) +/2)i;. Unlike 
the earlier example it is difficult to integrate this equation further and obtain the 
general solution. Therefore one can look for variables which transform the third 
order nonlinear ODE (5.7) to a third order linear ODE so that the non-trivial 
integration can be avoided. Now using the equation (4.6) we get w = J ^dt = t—j. 
The new variables 

w= — : — , 2; = logo;, (5.9) 

X 

transform the equation (5.7) to the form (4.9). Unlike the earlier example, w and z 
admit the variable x explicitly and so they become contact transformation for the 
given equation. 

Integrating the hnear third order equation (4.9), we get w = ^z^ + I2Z + I3, 
where Ii, i = 1,2,3, are integration constants. Now replacing w and z in terms of 
the old variables and using the previous result x = (/i log(a;) + /2)a;, one can obtain 
the general solution for the equation (5.7) in the form 



x{t) = ( - Ji ± J If + Ii - 2h{h - t) ) e ^ . (5.10) 
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(c) Example 3: Generalised Sundman transformation 

Next we consider the hydro dynamic type equation of the form (Berkovich & 
Orlova 2000; Euler & Euler 2004) 

X = — , (5.11) 

X 

which admits a first integral in the form 7i = | and the latter can be rewritten 
as I\ = = from which wc identify w = x and dz = xdt. By utilizing the 

new variables, one can transform (5.11) to the second order free particle equation, 
— 0. However, from equation (4.6) we get w = J xxdt = x^. Then the Sundman 
transformation 



X 



^ dz = xdt, (5.12) 



transforms the equation (5.11) to the form (4.9), namely = 0. 

To derive the solution we proceed as follows. Rewriting the first integral Ii in 



X dt ' 

in the latter expression we get w = hz. From equation (4.6) (for the present case 

dz ^ 
dt ' 



the integral form we get, h = ^ x = Ii J xdt. Now using the identity (5.12) 

on we get 

Gi = w and G2 — %) we have 



w ■ 



J wdz = J hzdz = ^z^ + I2, (5.13) 



where I2 is the integration constant. Using (5.12) in (5.13) we obtain x"^ = ^z^+l2, 
which in turn leads to a diff'erential equation which connects the variables z and t 
in the form (using the relation dz = xdt) 



dz = ^ ^z^ + hdt. (5.14) 

Integrating (5.14), we obtain z = y^(e^(*+^3) _ g-Vhit+h)^^ where I3 is the 

integration constant. Substituting the latter in the relation = ^z"^ + I2, we 
arrive at the general solution for (5.11) in the form 

x{t) = ^/hcOsh^/h{t + I3). (5.15) 

We note that the solution for the equation (5.11) has been already derived in an 
alternate way from the Sundman transformation (5.12) by Euler & Euler (2004). 
However the procedure which we described in the above is new and can be used for 
more general linearizing transformations also, as we see below. 

(d) Example 4: Generalized linearizing transformation (GLT) 

As we noted earlier some nonlinear ODEs can be linearized only through more 
general nonlocal form of transformations, which we designate here as generalized 
linearizing transformations (GLTs). We illustrate the GLT with the same example 
discussed as Example 2 in §3, which admits a first integral of the form Ii = ^ (vide 



Article submitted to Royal Society 



Integrability and Linearization 



19 



equation (3.16)). Rewriting this first integral in the form (4.2) we get Ii = fj^(a;), 
so that 

X 

w = X, dz = —dt, (5.16) 

X 

which can be effectively used to transform the nonlinear ODE (3.16) to the equation 
= 0. Using the equation (4.6) we get w = J xdt = x. Then the GLT becomes 

w = X, dz = —dt, (5.17) 

X 

which can be used to transform the equation (3.16) to the form w"' = 0. Note that in 

the present case the new independent and dependent variables admit x and x terms, 
respectively and that the transformation is nonlocal. Indeed no such linearizing 
transformations have been reported in the literature atleast to our knowledge. We 
also now establish a method of finding the general solution for this case. 

,2 

Integrating once the equation = we get w = Iiz from which we obtain 

X = hz. (5.18) 

On the other hand equation (4.6) provides us a relation (after using (5.16) and 
(5.17)) 

x = ^^z^ + l2. (5.19) 



Now using (5.18) in (5.19), we obtain 

2/i 



dz = dt. (5.20) 



JlZ^ + 2l2, 

The variables are now separated out and one can integrate (5.20) and obtain 

^= \/¥^tanW:^(< + /3). (5.21) 
V Ji \ ^1 

Now substutiting (5.21) into (5.18) we get 

x{t) = tan \h^it + h). (5.22) 

V -fi 

which is effectively the same as (3.25). 

Finally, we note that the procedure given above can be profitably utilized for 
other examples also which are linearized by GLTs. 



(e) Exam,ple 5: An elementary nontrivial system of hydrodynamic type 

Finally, to show the importance of the GLT and how this transformation gives 
additional information about the linearization of nonlinear third order ODEs we 
consider the following specific example which was discussed in the literature (Berkovich 
1996; Berkovich & Orlova 2000), 

'x = — — 4ax'^x, a : parameter. (5.23) 

X 
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Equation (5.23) is nothing but the dynamical equation of the Euler-Poinsot case 
of a rigid body written in terms of a single variable (Berkovich 1996; Berkovich & 
Orlova 2000) . As wc have seen earlier this equation is linearizablc in the case a = 
through generalized Sundmann transformation. However, we wish to show here that 
the general equation (5.23) itself is linearizable through the GLT. 

From the first integral /i = - + 2ax'^ associated with equation (5.23), one can 
identify the GLT 

w = x'^, dz = dt, (5.24) 

which transforms the equation (5.23) to the form (4.9). Note that for the choice a = 
the independent variable becomes dz = 2xdt and so it becomes the generalized 
Sundman transformation, equation (5.12), identified in the literature (Berkovich & 
Orlova 2000; Euler & Eulcr 2004). Now following the steps given in Example 4 one 
can deduce the general solution for the equation (5.23) in terms of Jacobian elliptic 
function as 

x{t)= (ji{c - (c - b)sn^yah{c - a){t - to),m]) + , (5.25) 

wherea + 6 + c= ^2£l_3£M, ^6 + ac + c6 = M^^, a6c = = ^ 
and I2 and to are integration constants. 



6. Conclusion 

In this paper we have discussed a method of finding the integrals of motion and 
general solution associated with third order nonlinear ODEs through the modified 
PS method by a nontrivial extension of our earlier work on second order ODEs 
(Chandrasekar et al. 2005). We illustrated the validity of the method with suit- 
able examples. Further, we introduced a technique which can be utilized to derive 
linearizing transformations from the first integral. Interestingly, wc showed that dif- 
ferent types of transformations, namely, point, contact, Sundman and generalized 
linearizing transformations can be derived in a unique way from the first integral. 
We also indicated a procedure to derive general solution for the third order ODEs 
when GLTs occur. We believe that the GLT introduced in this paper will be highly 
useful to tackle new systems such as equation (5.23). Finally, the modified PS 
method can also be extended to higher order ODEs and coupled systems of ODEs. 
As far as the linearization of higher order ODEs are concerned it is still an open 
and challenging problem. As we pointed out in the introduction one can unearth 
a wide variety of linearizing transformations for the higher order ODEs besides 
formulating the necessary and sufficient condition for linearizing these equations in 
each form of transformation. We hope to address some of these aspects shortly. 

The work of VKC is siipportcd by Council of Scientific and Industrial Research 
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